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Section  1 


1.0  Executive  Summary 

The  objective  of  this  work  is  the  development  of  a  process  to  fabricate  an  axial  grooved 
heat  pipe  embedded  in  printed  wiring  board  (PWB)  technology  and  to  evaluate  its 
thermal,  mechanical  performance,  and  reliability.  The  results  clearly  show  that  the 
embedded  heat  pipe  was  operable  and  the  heat  transfer  through  the  substrate  was 
enhanced  (see  section  3).  However,  the  heat  pipe  was  limited  in  the  power  input  it  could 
handle,  as  heat  pipe  failure  occurred  around  10  W  by  delamination  of  the  PWB  and 
subsequent  loss  of  hermeticity  and  loss  of  the  working  fluid.  The  failure  was  detected 
visually.  After  failure  occurs,  the  heat  pipe  begins  to  behave  as  an  uncharged  heat  pipe. 
The  temperatures  in  the  whole  PWB  starts  to  increase  until  it  reaches  the  temperature 
distribution  of  the  uncharged  heat  pipe  for  the  corresponding  input  power. 


Temperature  behavior  after  failure  at  point  Tl. 
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The  development  of  a  hermetic  structure  within  the  heat  pipe  structure  to  contain  the 
working  fluid  required  development  of  a  plating  system  that  allows  plating  within  the 
channel  structure.  Since  the  operation  of  the  wick  in  the  heat  pipe  is  strongly  dependent 
on  the  ability  of  the  working  fluid  to  wet  the  wick,  the  wetting  angle  of  the  plating 
material  was  measured  for  the  various  materials  options  available  (bare  laminate,  copper 
plating,  other  plating).  After  the  copper  plating  the  wetting  angle  of  the  inner  surface  of 
the  embedded  heat  pipe  has  been  reduced  from  104  degrees  (FR404  Epoxy  with  DI 
water)  to  56  degrees  (plated  copper  with  DI  water). 

Plating  chemistries  are  available  from  the  PWB  industry  including  both  electroless  and 
electroplating  solutions.  To  accomplish  this,  a  microplating  system  was  fabricated  which 
pumped  plating  solutions  thru  the  microchannel  from  a  recirculating  chamber.  The 
microplating  container  was  designed,  with  control  on  flow  speed  and  anode  position. 
Additionally,  the  pumping  mechanism  was  chosen  to  ensure  the  fluid  in  the  microchannel 
is  not  electrically  isolated  from  the  plating  container.  Five  configurations  of  plating  cells/ 
structures  were  evaluated. 

The  mechanical  properties  of  the  copper  plated  film  was  determined  by  tensile  testing  as 
these  properties  were  required  to  model  the  mechanical  performance  (i.e.  thermally 
induced  stresses  and  strains)  over  the  environmental  stress  exposure.  The  approach  to 
measuring  the  stress-strain  curve  is  similar  to  the  standardized  ASTM  E  345-86.  The 
electrodeposited  copper,  which  is  4.7  mils  thick,  has  brittle  material  behavior.  The 
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copper  tested  here  shows  that  the  initial  elastic  portion  of  the  stress- strain  curve  is  not 


linear  and  to  determine  the  modulus  of  elasticity  the  secant  modulus  method  was  used. 

The  modulus  of  elasticity  measured  perpendicular  to  the  deposition  direction  is  5720  ksi, 
which  is  lower  than  the  typical  material  value  (16700  ksi).  It  has  a  large  standard 
deviation  (13  %). 

The  yield  strength  (12.946  ksi)  and  percent  of  elongation  (48)  are  very  close  to  the  typical 
material  values  (10  ksi  and  45  respectively).  The  typical  material  value  is  a  factor  of  two 
larger  than  the  fracture  strength. 


Comparison  of  results  from  tensile  test  and  typical  properties 


E 

(psi) 

sy 

(psi) 

Sf 

(psi) 

Elongation 

(%) 

Experimental  Properties 

5.720  x  106 

12946 

16775 

48 

Typical  Properties 

16.7  x  106 

10000 

34075 

45 

The  results  of  tension  tests  from  selected  portions  of  the  electrodeposited  copper  material 
may  not  totally  represent  the  strength  and  ductility  of  the  entire  plated  surface  of  the 
embedded  heat  pipe  of  its  in-service  behavior  in  different  environment. 

The  physical  structure  of  the  embedded  heat  pipes  was  analyzed  by  finite  element 
analysis  (FEA).  (  see  section  4)  Figure  4.  12  to  Figure  4.17  compare  the  FEM  results 
with  the  experimental  data.  The  thermal- stress  finite  element  analysis  of  the  embedded 
heat  pipe  resulted  in  data  that  provide  information  of  the  strain  and  stress  generated  in  the 
substrate  for  different  level  of  input  power  at  different  points  along  the  heat  pipe,  with 
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this  information  the  temperature,  stress,  and  location  were  the  heat  pipe  could  fail  is 
inferred. 

1.1  Overview  of  Micro  Heat  Pipe  Technology 

Thermal  management  is  becoming  a  limiting  factor  in  system  applications  due  to  the 
increased  power  of  semiconductor  device.  The  Semiconductor  Industries  Association 
(SI A)  roadmap  projects  operating  power  level  of  advanced  microprocessors  at  over  200 
W/die  by  2005.  This  is  further  compounded  by  the  increased  complexity  at  the  system 
level  and  the  development  of  new  high  power  technologies,  such  as  SiC,  which  are 
projected  to  operate  at  power  levels  far  exceeding  current  devices.  These  requirements 
are  placing  thermal  demands  that  cannot  be  handled  by  current  technology  or  require 
technologies  that  add  great  cost  and  weight  at  the  system  level. 

To  provide  thermal  management  for  high  heat  dissipating  devices  mounted  on  laminate 
boards  when  air-cooling  is  not  an  option  requires  that  a  conduction  path  for  the  dissipated 
heat  must  be  provided  directly  from  the  device  to  the  sink.  For  power  devices,  this  is 
typically  accomplished  by  bonding  a  metal  heat  sink  to  the  circuit  card  and  mounting  the 
component  directly  on  the  heat  sink.  For  the  heat  sinks  to  work  adequately,  the  heat  sinks 
are  usually  3  mm  thick  or  greater,  and  need  to  be  machined  and  finished  and  the  heat  path 
still  must  travel  down  the  length  of  the  heat  pipe  to  the  thermal  sink  unless  direct  fluid 
cooling  is  used.  With  the  project  heat  fluxes  of  new  high  power  components,  conduction 
thermal  paths  will  not  be  able  to  handle  the  thermal  loads  without  incurring  a  large 
temperature  rise  from  the  component  to  the  heat  sink. 
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Previous  works  in  Florida  International  University  has  addressed  the  need  for  system 
level  integration  of  thermal  management  tools  by  the  development  of  micro  heat  pipes 
fabricated  within  a  ceramic  co-fired  interconnect  substrate.  Micro  heat  pipes,  capable  of 
transporting  large  amounts  of  heat  (40  W)  over  relative  long  distances  (10  cm 
demonstrated)  have  been  fabricated  in  both  high  temperature  alumina  and  low 
temperature  cofire  ceramic  (HTCC  and  LTCC).  The  heat  pipe  was  found  to  have  an 
effective  thermal  conductivity  greater  than  10,000  W/m-°K,  approximately  five  times 
than  diamond  and  over  1000  times  higher  than  the  ceramic  it  replaced  within  the 
substrate.  This  allows  an  entirely  new  set  of  thermal  management  options  in  a 
technology,  such  as  LTCC,  where  the  thermal  conductivity  of  the  material  is  less  than  4 
W/m-°K. 

Heat  pipes  are  designed  to  transport  heat  over  relatively  long  distances  with  no  additional 
input  power  (i.e.  passive  function).  A  heat  pipe  consists  of  a  sealed,  hermetic  enclosure, 
with  three  distinct  regions:  an  evaporator,  a  condenser,  and  an  adiabatic  region  separating 
these  two  regions.  The  enclosure  contains  a  working  fluid,  which  absorbs  heat  by 
evaporation  at  the  evaporator,  travels  as  a  vapor  in  the  adiabatic  region  to  the  condenser, 
where  the  heat  is  removed.  The  working  fluid  returns  to  the  evaporator  section  by  the 
capillary  action  of  a  wick  structure.  The  total  heat  energy  (W  or  J/s)  that  the  pipe  could 
transport  relates  to  the  heat  of  vaporization  (J/kg)  and  the  amount  of  fluid  undergoing 
evaporation  (kg/s).  The  heat  pipe  used  an  axially  grooved  wick  structure  to  provide 
transport  of  the  working  fluid  from  the  condenser  to  the  evaporator.  The  design  of  the 
wick  is  the  most  critical  aspect  to  insure  high  fluid  transport  and  thus  high  heat  transport 
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capability.  Different  designs  of  wicks  have  been  developed  by  researchers,  the  most 


common  types  of  homogeneous  wicks  included  wrapped  screen,  sintered  metal,  and 


small  axial  grooves  along  the  inner  walls  of  the  heat  pipe  [1]. 
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Figure  1.1  Heat  pipe  structure  and  operation. 


Initial  efforts  in  the  oeramic  heat  pipes  used  water  as  the  working  fluid,  although  other 
working  fluids  are  also  available  based  on  the  desired  operating  temperature  and  other 
systems  requirements  [2], 

The  most  significant  limiting  factor  in  the  power  handling  capability  of  miniature  heat 
pipes  is  the  capillary  limit,  which  is  the  point  at  which  the  wick  cannot  deliver  working 
fluid  at  the  same  rate  at  which  it  is  evaporated.  This  condition  leads  to  evaporator  dry¬ 
out,  indicated  by  a  rapid  increase  in  the  local  temperature  of  the  evaporator.  The 
capillary  limit  is  a  function  of  the  surface  tension  of  the  working  fluid  and  the  frictional 
forces  caused  by  the  interaction  of  the  working  fluid  with  the  wick  material  and  with  the 
opposing  vapor  flow.  For  axially  grooved  heat  pipes,  the  capillary  limitation  has  been 
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well  developed  by  others  [3],  and  has  been  reviewed  for  miniature  embedded  heat  pipes 
by  [4], 

The  design  of  an  optimized  miniature  heat  pipe  for  the  thermal  management  of 
electronics  involves  a  complex  balance  of  system  requirements,  and  the  thermophysics  of 
heat  pipes,  which  involve  capillary  forces,  phase  change,  two -phase  flow  interaction,  and 
porous  media  and/or  channel,  flow.  The  relevant  issues  will  be  discussed  starting  from 
the  evaporator  end  of  the  heat  pipe  where  the  heat  is  generated,  and  will  end  at  the 
condenser  end  where  heat  is  removed. 

1.1.1  The  Capillary  Limit 

To  understand  the  capillary  limitation  with  respect  to  both  width  (Wg)  and  depth  (Dg)  of 
straight  rectangular  grooves  for  a  given  internal  height  and  width  of  the  housing,  it  is 
opportune  to  refer  to  the  following  analysis  developed  by  Zampino  [2], 


Figure  1 .2  Heat  pipe  cross  section. 


7 


The  capillary  limit  is  a  function  of  the  surface  tension  of  the  working  fluid  and  the 
frictional  forces  caused  by  the  interaction  of  the  working  with  the  wick  material  and  with 
the  opposing  vapor  flow.  For  axially  grooved  heat  pipes,  the  capillary  limitation  is  given 
in  Equation  1.1  [3],  which  relates  the  pumping  action  of  the  working  fluid  due  to  surface 
tension,  a,  and  the  frictional  forces,  F,  and  Fv ,  to  flow  of  the  working  fluid  from  the 

condenser  back  to  evaporator.  The  frictional  forces  are  described  by  two  components,  Ft , 
which  represents  the  interaction  of  the  working  fluid  and  the  wick  material  (structure), 
and  Fv ,  which  represents  the  resistance  to  the  vapor  flow  due  to  the  geometry  of  the 
vapor  space.  Leff  is  the  effective  length  of  the  heat  pipe,  which  is  equal  to  La  +  0.5  (Le  + 
Lc). 


Qcap  ~ 


2 c 


WgLeff(Fl+Fv) 


(1.1) 


The  liquid  and  vapor  frictional  coefficient  are  given  in  Equations  (1.2)  and  (1.3). 
Referring  to  Equation  (1.2),  it  is  seen  that  the  permeability  of  the  wick,  Kg ,  is  inversely 

proportional  to  the  liquid  frictional  coefficient.  Aw  is  the  cross-sectional  area  of  the  wick 
structure.  Referring  to  Equation  (1.3),  it  is  seen  that  the  vapor  frictional  coefficient  is 
proportional  to  the  friction  factor,  /(Re vh)  which  is  given  in  equation  (1.7)  [5], 


F,= 


P, 


PiAvKghfg 


(1.2) 
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The  hydraulic  diameter  for  the  vapor  flow  space  Dv>h=2Rv>h  ,  and  that  of  the  axially 
grooved  wick  structure  D^h  ,  can  be  evaluated,  respectively,  by  the  following  relations: 

Dv,h  =  2  tvWv/(tv  +  Wv)  (1.4) 

Di,h  =  4DgWg/(2Dg  +  Wg)  (1.5) 

The  permeability  of  the  axial  grooved  wick  is  given  in  Equation  (1.6),  which  shows  a 
ratio  between  the  flow  area  and  the  friction  to  the  flow.  The  friction  factor,  (/ Re//r) ,  is 

determined  using  Equation  (1.8)  and  (1.9)  which  include  the  shear  stress  interaction 
between  the  vapor  flow  and  the  liquid  [6], 


K  _  Dj$_ 

S  2(/Rea) 


(1.6) 


where  <|)  can  be  calculated  by  <j)  =  Wg  /  Sg 

(/Rev h)  =  24(1-  1.3553a  +1.9467a2  -1.7012a3  +0.9564a4  -0.2537a5)  (1.7) 

where  a  =  Wv/tv,  when  Wv  <  tv  or  a  =  tv/Wv ,  when  Wv  >  tv 
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All  the  equations  thus  far  show  the  paradox  of  miniature  heat  pipe  design,  that  is,  the  heat 
transport  capability  is  increased  by  increasing  the  flow  rates  of  vapor  and  fluid,  however, 
large  flow  rates  through  small  channels  leads  to  high  frictional  losses  which  decrease  the 
heat  transport.  Hence,  the  optimization  of  the  design  parameters  must  be  performed  to 
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provide  an  efficient  design.  Insight  into  an  optimal  design  can  be  gained  by  showing  the 
effect  of  the  groove  depth,  Dg,  and  width,  Wg  (see  Figure  1.6),  on  the  capillary  limit  as 
shown  in  Figure  1.3.  From  this  figure  it  is  clear  that  an  optimal  groove  width  exists 
which  is  approximately  the  same  for  various  groove  depths  (all  other  parameters  held 
constant),  and  as  the  groove  depth  increases,  so  does  the  capillary  limit. 


-*-A)  Dg  =  0.64  mm 
-“-B)  Dg  =  0.38  mm 
C)  Dg  =  0.25  mm 


Figure  1.  3  The  effect  on  the  capillary  limit  due  to  variation  in  groove  width  for  groove 
depths  of:  (A)  0.64  mm,  (B)  0.38  mm,  and  (C)  0.25  mm. 


In  an  effort  to  minimize  the  thickness  of  the  heat  pipe,  it  is  desirable  to  minimize  the 
height  of  the  vapor  space.  The  effect  of  the  vapor  space  height  is  shown  in  Figure  1 .  4, 
which  shows  that  the  capillary  limit  initially  increases  sharply  with  increasing  vapor 
space  height  and  then  becomes  relatively  constant.  Again,  the  effect  of  the  groove  depth 
is  seen  and  follows  the  same  trend  as  discussed  for  Figure  1.3.  An  important  conclusion 
drawn  from  the  figure  is  that  the  minimum  vapor  space  height  is  closely  related  with  the 
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depth  of  the  grooves,  and  that  additional  vapor  space  height  will  not  have  significant 


effect. 


A)  Dg  =  0.64  mm 
-°-B)  Dg  =  0.38  mm 
C)  Dg  =  0.25  mm 


Figure  1.  4  The  effect  on  the  capillary  limit  due  to  variations  in  vapor  space  height  for 
groove  depths  of:  (A)  0.64  mm,  (B)  0.38  mm,  and  (C)  0.25  mm. 

The  effect  of  the  interaction  between  the  groove  width  and  the  vapor  space  height  is 
made  clearer  by  holding  the  depth  of  the  grooves  constant  and  varying  the  other  two 
parameter  as  shown  in  Figure  1.5.  The  curves  in  the  figure  show  that  for  small  vapor 
space  heights,  smaller  groove  widths  provide  the  highest  capillary  limit,  but  as  the  height 
of  the  vapor  space  increases,  increases  in  the  capillary  limit  are  obtained  with  larger 
groove  widths.  The  “crossing- over”  effect  shown  by  the  curves  indicates  that  there  is  a 
competition  between  the  shear  stress  interaction  between  the  vapor  and  liquid  flows  and 
the  ability  of  the  wick  to  provide  adequate  mass  of  working  fluid. 
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tv  (mm) 


-*-  A)  Wg  =  0. 10  mm 
—  B)  Wg  =  0.25  mm 
C)  Wg  =  0.38  mm 


Figure  1.  5  The  effect  on  the  capillary  limit  due  to  variations  in  vapor  space  height  for 
groove  widths  of:  (A)  0.10  mm,  (B)  0.25  mm,  and  (C)  0.38  mm. 

Small  vapor  space  heights  cause  greater  shear  stress  interaction.  Hence,  only  smaller 
groove  widths  can  minimize  the  surface  area  of  the  working  fluid  interacting  with  the 
vapor  flow.  However,  at  some  point,  when  vapor  height  exceed  certain  value  the  shear 
stress  interaction  will  become  negligible  and  the  capillary  limit  becomes  strongly 
dependent  on  the  cross-sectional  area  through  which  the  working  fluid  can  flow. 

1.1.2  Surface  Tension  and  Surface  Wetting 

As  discussed  earlier,  the  operation  of  the  wick  in  the  heat  pipe  is  strongly  dependent  on 
the  ability  of  the  working  fluid  to  wet  the  wick.  The  better  the  working  fluid  wets  the 
wick  material,  the  greater  the  surface  tension  forces  will  be  on  the  working  fluid.  The 
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capillary  limit  as  described  is  a  balance  between  the  surface  forces  due  to  capillary  action 
and  the  retarding  forces  due  to  gravity  and  friction.  Therefore,  any  enhancement  of  the 
wick’s  surface  that  increases  the  wetting  of  the  working  fluid  will  increase  the 
performance  of  the  heat  pipe. 

The  most  basic  model  for  determining  the  surface  tension  forces  was  introduced  early  in 
this  century  by  Washburn  and  Rideal  [7,  8]  as 

Fa  =  27tra  cos(0 )  (1.10) 

This  equation  shows  the  relationship  of  the  main  parameter  affecting  the  wetting  of  a 
fluid  to  a  capillary  groove,  specifically,  the  size  of  the  capillary,  the  surface  tension,  and 
the  wetting  angle.  The  capillary  radius  is  mainly  driven  by  the  geometry  of  the  heat  pipe 
and  the  fabrication  technologies  used.  Likewise,  the  selection  of  working  fluid  is  limited 
to  a  few  fluids  that  are  compatible  with  the  environmental  requirements  of  the  system. 
Therefore,  the  only  parameter  left  which  can  be  optimized  is  the  wetting  angle,  which  is  a 
function  of  the  working  fluid  and  the  heat  pipe  material  selection.  Again  the  selection  of 
the  heat  pipe’s  material  depends  on  the  fabrication  technologies  employed  and  is  limited. 
However,  surface  treatments  can  be  employed  to  the  selected  materials  to  affect  the 
wetting  angle.  Liquids  such  as  water  have  relatively  low  wetting  angles  on  materials 
such  as  gold  (approximately  12  degrees).  From  equation  1.10  it  can  be  seen  that  small 
wetting  angles  will  increase  the  surface  tension  force. 

An  advantage  that  laminate  technology  has  over  ceramic  is  the  wider  choice  of  materials 
that  can  be  plated.  Ceramic  technology,  due  to  the  need  for  high  temperature  firing,  is 
limited  to  metal  organics,  which  can  be  used  to  coat  the  channel.  However,  plating 
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compositions  compatible  with  printing  wiring  board  (PWB)  technology  include  the 
precious  metals  (Au,  Pt,  Ag,  Pd,  Rh,  etc.),  metals  such  as  Ni,  Cr,  Pb,  Sn,  Zn,  etc.,  and 
their  alloys.  This  provides  a  wide  range  of  possible  materials  to  optimize  wetting 
performance  in  this  technology. 

As  mentioned  previously,  the  predominate  technology  for  interconnect  boards  is  laminate 
technology.  This  work  considers  extending  the  development  of  micro  heat  pipes  into 
laminate  interconnect  board  (printed  wiring  board).  A  major  drawback  of  laminate 
technology  is  the  inherent  low  thermal  conductivity  of  the  polymer  laminate  materials, 
typically  in  the  range  of  0.1  W/m°K.  This  development  will  utilize  the  laminate  structure 
of  printed  wiring  boards  to  develop  the  wick  structure. 

The  vapor  space  (tv)  for  the  heat  pipe  can  be  easily  formed  in  the  laminate  stack  by 
punching  slots  in  each  layer  of  material.  Stacking  several  layers  of  laminate  material 
with  the  slots  cut  into  them  forms  a  long  rectangular  duct  within  the  substrate.  The 
fabrication  of  the  axial  grooves  for  the  wick  structure  is  shown  schematically  in  Figure 
1.6. 


14 


Figure  1.6  Schematic  representation  of  heat  pipe  section  showing  axial  groove 
construction 

Axial  grooves  in  the  top  and  bottom  interior  surface  of  the  heat  pipe  could  be  fabricated 
by  numerically  machining  small  grooves  into  the  laminate  layers  of  the  top  and  bottom 
layers  of  the  microchannel  or  plating  up  grooves  in  these  copper  coated  laminate  layers. 
The  first  process  requires  the  use  of  micro  end-mills  and  numerically  controlled  three 
axis  machining  system  with  a  precision  of  at  least  0.025  mm,  while  the  second  requires 
relatively  thick  plating  structures  with  fine  line  resolution.  In  contrast,  the  fabrications  of 
axial  grooves  in  the  sidewalls  of  the  heat  pipe  are  relatively  easy.  Varying  the  width  of 
the  rectangular  slot  cut  into  each  laminate  layer  to  form  he  vapor  space  forms  these 
grooves.  The  variation  in  the  width  of  the  slots  on  alternating  layers  provides  control 
over  the  depth  ( Dg )  of  grooves.  However,  the  width  (WH)  of  each  groove  is  dictated  by 
the  thickness  of  the  tape  layer.  Fortunately,  typical  laminated  material  is  0.1  to  0.5  mm 
(0.004  to  0.020  in)  in  thickness,  which  allows  for  fine  narrow  grooves  to  be  formed. 
Both  types  of  grooves  have  been  demonstrated  in  laminate  construction  in  alumina  co- 
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fire  ceramic  materials.  However,  the  ease  of  manufacturing  with  sidewall  grooves 
greatly  reduces  manufacturing  cost  with  minimal  performance  impact. 

The  heat  pipe  will  be  assembled  using  patterned  b- stage  prepreg,  which  will  be  laminated 
using  standard  printed  wiring  board  manufacturing  processes.  Initially,  it  is  anticipated 
that  the  heat  pipe  will  be  manufactured  at  the  completion  of  the  PWB  manufacturing 
(thru  hole  plating,  solder  mask,  solder  plating,  OSP  deposition,  etc).  It  is  envisioned  that 
capture  pads  for  the  fill/scaling  arrangement  would  be  post  drilled  at  both  ends  and  the 
microchannel  subsequently  plated. 

Section  2 

Micro  Heat  Pipe  Fabrication 

The  heat  pipe  is  assembled  using  FR404  epoxy  laminate  layers  and  B15  prepreg  bonding 
sheets,  which  is  laminated  using  standard  printed  wiring  board  manufacturing  processes. 
Initially,  it  is  anticipated  that  the  heat  pipe  will  be  manufactured  at  the  completion  of  the 
PWB  manufacturing  (thru  hole  plating,  solder  mask,  solder  plating,  OSP  deposition,  etc.). 
It  is  envisioned  that  capture  pads  for  the  charging/sealing  arrangement  would  be  post 
drilled  at  both  ends  and  the  micro  channel  subsequently  plated. 

Figure  2.1  shows  the  cross  sectional  view  of  the  side  walled  grooved  laminate  heat  pipe. 
The  side  waned  grooves  for  the  heat  pipe  are  formed  in  the  laminate  stack  by  punching 
slots  in  each  laminate  layer  of  FR404  material.  Stacking  several  layers  of  laminate 
material  with  the  slots  cut  into  them  forms  a  long  rectangular  duct  within  the  laminate 
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substrate.  Varying  the  width  of  the  rectangular  slot  cut  into  each  layer  to  form  the  vapor 
space  for  these  grooves.  Variation  in  the  width  of  the  slots  on  alternating  layers  provides 
control  over  the  depth  (Dg)  of  the  grooves.  The  width  ( Wg)  of  the  grooves  is  controlled  by 
the  thickness  of  the  tape  layers  (see  Figure  2.1).  The  FR404  epoxy  laminate  layers  with 
the  slots  cut  into  them  were  bonded  together  during  lamination  process  by  B15  prepreg 
bonding  sheet  to  form  the  heat  pipe. 

Fabrication  of  the  heat  pipe  involves  critical  steps  including  etching,  lamination, 
electroplating,  and  fluid  charging. 


TOP  LAYER 
BONDING  SHEET 
LARGER  SLOT  LAYER 


SMALLER  SLOT  LATER 


BDTTOM  LAYER 


|— —  Dg 


Wg 


Figure  2.1  Laminate  layers  setup. 

The  laminate  layers  itself  and  its  rectangular  slots  were  cut  using  a  numerically  controlled 
three-axis  machine  system  (OZO  system).  The  OZO  system  supports  the  ANSI/IPC-NC- 
349  computer  numerical  control  format  for  drills  and  routers.  See  appendix  B  for  the 
cutting  computer  program.  It  is  necessary  that  original  prepreg  bonding  sheets  have  a 
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larger  slot  than  the  larger  slot  layer  (see  Figure  2.2)  because  during  the  lamination  the 
bonding  sheet  material  will  flow  and  reach  its  final  dimension  similar  to  the  larger  slot 
layer  dimension  (see  Figure  2.1). 

It  is  important  to  state  that  US  Units  are  used  in  the  rest  of  this  work  but  keeping  watts 
(W)  as  the  power  unit,  since  it  is  the  power  unit  most  generally  used  in  the 
microelectronics  industry. 


Figure  2.2  Laminate  layers  dimensions  (inches):  (a)  smaller  slot  layer  (17  mils  thick), 
(b)  larger  slot  layer  (7  mils  thick),  (c)  bonding  sheet  (3  mils  thick). 


2.1  Etching 

Etching  of  the  FR404  laminate  layers  is  necessary  to  ensure  good  bonding  to  the  prepreg 
sheets  during  lamination.  The  laminate  layers  to  be  etched  are  immersed  in  copper 
cleaner  to  remove  any  surface  contaminants  like  dust,  stains,  organic  particles,  etc.,  that 
may  have  deposited  during  handling  and  routing.  These  surface  contaminants  result  in 
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poor  adhesion  of  the  photoresist  during  the  photolithography  process  and/or  affect  the 
subsequent  processes  like  etching  and  lamination  [9]. 


Figure  2.3  FR404  Laminate  layer  after  etching. 

The  master  artwork  of  the  desired  pattern  is  generated  using  computer-aided  design 
(CAD)  package.  The  artwork  generated  was  transferred  on  to  a  transparency,  which 
served  as  the  mask. 

The  laminate  layers  are  patterned  using  standard  photolithography  process  for  selective 
removal  of  materials.  The  Patterning  is  done  by  spinning  a  film  of  KPR-3  photoresist  on 
the  surface  of  FR404  laminated  layers.  The  speeds  ranged  from  3000  to  4000  rpm,  40 
sec.  An  alternate  method  to  obtain  a  uniform  coat  of  photoresist  on  the  part  is  by  dipping 
the  part  in  the  photoresist  and  drawing  it  out  of  at  a  constant  rate  (approximately  lin/min) 
[10]. 
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The  resist- coated  part  is  soft  baked  for  5  min.  at  250  °F  in  a  conventional  oven.  The 
purpose  of  the  soft  bake  step  is  to  removes  any  solvents  on  the  surface  and  dry  the 
photoresist. 

The  negative  artwork  and  a  negative  photoresist  were  used  for  the  fabrication  process. 
Negative  acting  photoresists  polymerize  when  exposed  to  an  ultra  violet  light  source. 
The  artwork  is  sandwiched  between  the  piece  of  glass  and  copper  clad  board  and  exposed 
to  an  ultra-violet  light  source  through  a  photo  mask.  For  good  results  the  exposure  time 
lies  between  4  and  5  min. 

The  parts  are  dipped  in  a  developer  to  get  a  precise  definition  of  the  desired  image.  Then 
laminated  layers  are  rinsed,  dried  and  post-  baked  at  250  °F  for  10  min. 

The  patterned  FR404  laminated  layers  are  etched  at  room  temperature  (70°F)  using 
chemical  etching.  The  patterned  board  is  immersed  in  the  Ferric  Chloride  solution 
(copper  etchant)  for  25  to  40  min.  Continual  agitation  is  essential  to  achieve  a  uniform 
etching  rate.  The  board  is  immersed  in  stripper  solution  to  remove  the  residual  photoresist 
film,  rinsed  in  cold  water,  and  dried  [10]. 

2.2  Lamination 

The  manufacturing  recommended  press  cycle  for  FR404  is  as  follows  [11]: 

1.  Pre-heat  plate  to  350  to  360  °F. 

2.  Load  and  center  all  layers  as  quickly  as  possible. 

3.  Close  press  to  300-350  psi. 

4.  When  center  temperature  of  staked  layers  reaches  350  °F,  cure  for  50  min. 
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5.  Cool  material  under  moderate  to  normal  pressure  slowly  (5  °F/min)  through  250  °F. 


Cooling  rate  can  be  increased  when  temperature  lower  than  250  °F. 

Since  the  press  time  is  depend  on  the  stack  height,  number  of  bonding  prepreg  sheets,  and 
temperature  used,  in  practical  application,  the  above  recommended  lamination  parameters 
can  be  varied  in  order  to  obtain  the  best  results  for  the  specific  application.  In  this  work, 
several  lamination  procedures  were  used  and  tested,  the  result  is  summarized  in  Table 
2.1. 


Table  2.1  Summary  of  the  lamination  process 


Lamination 

No. 

Temp. 

(°F) 

Time 

(min) 

Pressure 

(psi) 

Result 

1 

300 

60 

300 

Good  bonding,  bad  dimensional  accuracy 

2 

210 

180 

300 

Bad  bonding,  good  dimensional  accuracy 

3 

260 

60 

300 

Bad  bonding,  good  dimensional  accuracy 

4 

260 

180 

300 

Bonding  is  not  very  strong,  good 
dimensional  accuracy 

5 

280 

180 

300 

Bonding  is  not  very  strong,  good 
dimensional  accuracy 

6 

300 

120 

300 

Good  bonding,  bad  dimensional  accuracy 

7 

180 

300 

Good  bonding,  good  dimensional 
accuracy 
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Figure  2.4  PWB  with  embedded  heat  pipe. 


Figure  2.  5  External  dimensions  (inches)  of  the  PWB. 
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Figure  2.6  Dimensions  (inches)  of  the  PWB  cross-section  with  embedded  heat  pipe. 

2.3  Smear  Removal 

Routing  smear  refers  to  the  epoxy  resin  that  coats  the  inner  layer  copper  surface  in  the 
side  walled  grooves  formed  in  the  laminate  stack  and  occurs  during  the  routing  operation. 
This  smear  must  be  removed  before  plating  to  get  a  continuous  and  regular-coated 
metallization. 
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(a)  (b) 

Figure  2.7  Grooves:  (a)  without  smear  removal,  (b)  with  smear  removal. 


The  Sulfuric  Acid  method  was  used.  This  process  is  used  extensively  because  of  its  ease 
of  operation  and  reliability  of  results.  A  major  disadvantage  is  a  lack  of  control,  which 
leads  to  wall  pull  away.  Operator  safety  is  crucial  since  concentrated  sulfuric  must  be 
used.  The  method  process  is  as  follows  [12]: 

Sulfuric  Acid,  96  % 

20  s  dwell,  15  s  drain,  room  temperature 
Neutralizer 
8  min,  125°F 
Ammonium  bifluoride 
3  min,  room  temperature 
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High-pressure  hole  cleaning 


2.4  Copper  Plating 

To  make  the  microchannel  hermetic  and  have  enhanced  wetting  performance,  the  inner 
surface  of  the  channel  must  be  coated  with  a  metal  that  improves  these  properties.  This 
required  development  of  a  plating  system  that  allows  copper  plating  within  the  channel 
structure. 

Table  2.2  summarized  the  step  by  step  procedure  for  Conductive  Copper  Deposition 
processing  [13].  Each  step  is  named  and  has  its  respective  operational  time  listed.  As 
noted  in  the  chart,  a  total  of  61  minutes  is  required  to  go  through  this  process. 

The  experimental  testing  of  the  different  plating  setups  were  conduced  under  the 
following  conditions: 

Operating  temperature:  80  °F 

Solution  flow:  6  mL/min 

Current  density:  700  mA 

Agitation:  Solution  movement  (ultrasonic  vibration  bath) 


Tabb  2.2  Conductive  copper  deposition  process 


Step  No. 

Functions 

Time  Required 

Product  Name 

1 

Chemical  Clean  Holes 

5  Min. 

PTH  Cleaner 

2 

Water  Spray  Rinse 

1  Min. 

H20 

3 

Chemical  Clean  Copper 

30  Sec. 

Ammon.  Persol 

4 

Water  Spray  Rinse 

1  Min. 

h2o 
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5 

Acid  Tread  Panel 

5  Min. 

H2S04 

6 

Water  Spray  Rinse 

1  Min. 

H20 

7 

Acid  Tread  Panel 

5  Min. 

HCL 

8 

Activate- Seed  Hole 

10  Min 

Activator 

9 

Water  Spray  Rinse 

1  Min. 

h2o 

10 

Accelerate  Seeded  Palladium 

5  Min. 

Accelerator 

11 

Water  Spray  Rinse 

1  Min. 

H20 

12 

Overflow  Rinse 

1  Min. 

Deionized  H20 

13 

Electroless  Copper  Plate 

10  Min. 

Electroless  Cu 

14 

Water  Spray  Rinse 

1  Min. 

h2o 

15 

Acid  Treat  Panel 

1  Min. 

h2so 

16 

Electroplate  Copper 

10  Min 

Cu  Electrolyte 

17 

Water  Spray  Rinse 

1  Min. 

H20 

18 

Dry  Panel 

1  Min 

Air 

The  experimental  setup  for  different  systems  is  listed  as  follows: 

1.  Figure  2.8  shows  the  experimental  setup  for  electroplate  copper  in  a  closed  system 
with  the  anode  outside  of  the  heat  pipe 
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DC  Power  Source 


Figure  2.8  Experimental  setup  No.  1,  electroplate  copper  in  a  closed  system. 


Figure  2.9  Electroplate  setup  No.  1,  copper  plated  film. 

2.  Figure  2.  10  shows  the  experimental  setup  in  a  closed  system  with  the  anode  inside 
the  heat  pipe.  Prior  to  lamination,  the  copper  anode  is  inserted  inside  the  heat  pipe. 
Wax  is  used  as  an  insulator  to  separate  the  anode  with  the  cathode. 
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Figure  2.11  Electroplate  setup  No.  2,  copper  plated  film. 

3.  Figure  2.12  shows  the  experimental  setup  for  plating  in  a  closed  system  with  side 
holes  in  the  heat  pipe.  In  this  method,  the  copper  anode  is  inserted  into  the  heat  pipe 
through  a  hole  on  the  side  of  it,  and  some  insulating  glue  were  used  to  separate  the 
anode  from  the  heat  pipe.  The  electrolyte  is  pumped  through  the  heat  pipe  from  the 
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inlet.  After  plating,  the  side  holes  are  sealed  using  epoxy  materials,  although  metal 


plug  could  be  soldered  into  the  ends. 


Figure  2.12  Experimental  setup  No.  3,  electroplate  copper  in  a  closed  system  with  side 
holes. 


Figure  2.13  Electroplate  setup  No.  3,  copper  plated  film. 

4.  Figure  2.14  shows  the  experimental  setup  for  electroplate  copper  in  an  open  system. 
The  anode  is  a  rectangular  solid  copper  bar  inside  the  heat  pipe. 
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DC  Power  Source 


Figure  2.14  Experimental  setup  No.  4,  electroplate  copper  in  an  open  system. 


Figure  2.15  Electroplate  setup  No.  4,  copper  plated  film. 
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5  This  is  an  open  system  experimental  setup  (see  Figure  2.16).  The  only  difference 
between  this  setup  and  the  previous  setup  is  that  the  anode  is  a  rectangular  copper 
pipe  with  small  holes  on  it.  A  pump  is  used  to  pump  up  the  solution  through  the 
anode  to  inside  the  heat  pipe. 


Figure  2.16  Experimental  setup  No.  5,  electroplate  copper  in  an  open  system 


Figure  2.17  Electroplate  setup  No.  5,  copper  plated  film. 
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Section  3 


Processing  and  Testing 

3.1  Plated  Copper  Film  Tensile  Test 

The  method  used  covers  the  recording  of  the  complete  stress- strain  curve  of  copper 
plated  film  specimens  under  uniaxial  tension  at  room  temperature.  The  results  of  these 
tests  provided  the  information  on  modulus  of  elasticity  (E),  yield  tensile  strength  (sy), 
fracture  tensile  strength  (s  f),  and  elongation  (%EL)  of  the  copper  plated  film. 

3.1.1  Specimen  Preparation 

Two  different  size  specimens  were  fabricated  by  electroplating.  The  shape  of  the 
specimens  was  adopted  from  ASTM  E  345-86,  Standard  Method  of  Tension  Testing  of 
Metallic  Foil.  Five  samples  were  tested;  specimen  shape  No.  1  from  test  1  to  test  4  and 
specimen  shape  No.  2  in  test  No.  5.  The  dimensions  of  the  specimens  are  illustrated  in 
Figure  3.1. 
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Figure  3.1  Dimensions  of  the  copper  plated  film  specimens. 


3.1.2  Testing  Machine  and  procedure 

An  Instron  Mini-Tester  (1130)  with  500-Gram  Load  Cell  was  used  to  evaluate  the  stress- 
strain  behavior  of  the  plated  layer.  The  Mini- Tester  has  a  highly  sensitive  electronic  load 
weighing  system  with  load  cells  employing  strain  gages  for  detecting  the  load  applied  to 
the  specimen  under  test. 

The  specimen  is  fastened  between  two  grips  one  attached  to  the  load  cell  in  the  stationary 
crosshead  and  the  other  grip  attached  to  the  moving  crosshead  of  the  tester.  The  moving 
crosshead  is  operated  by  a  vertical  lead  screw.  A  synchronous  motor  drives  the  lead 
screw  through  a  gearbox  that  provides  a  variety  of  cross  head  speeds.  The  test  results  are 
present  on  a  strip  chart  recorder.  The  load  applied  to  the  specimen  under  test  is  recorded 
by  the  trace  of  a  pen  on  the  chart.  The  crosshead  displacement  is  recorded  by  the 
synchronous  motion  of  the  chart.  Thus,  the  chart  provides  an  accurate  recording  of  load 
versus  displacement. 
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3.1.3  Method  of  Determining  Elongation 


Chart  speed:  3  in./min 
Crosshead  speed:  0.01  in/min 
Gage  length:  1  in. 


w  ...  .  Chart  speed 

Magnification  = - 

Crossheat  speed 


(3.1) 


Chart  displacement  *  100 

hlongation{% )  = - 

Magnification  ratio* Gage  length 


(3.2) 


Table  3.1  Tensile  test  chart  displacements. 


Test  1 

Test  2 

Test  3 

Test  4 

Test  5 

Chart  Displacement 
(in) 

1.55 

0.925 

1.5 

1.35 

1.8 

Magnification 

300 

300 

300 

300 

300 

Elongation  (%) 

0.517 

0.3083 

0.5 

0.45 

0.6 

3.1.4  Mechanical  Properties 

Figure  3.2  shows  that  the  initial  elastic  portion  of  the  stress-strain  curves  are  not  linear. 
Hence,  it  is  not  possible  to  determine  a  modulus  of  elasticity  from  the  slope  of  the  linear 
segment.  For  this  nonlinear  behavior,  either  tangent  or  secant  modulus  is  normally  used 
[14]. 

In  this  analysis  the  secant  modulus  method  is  used.  The  secant  modulus  represents  the 
slope  of  a  secant  drawn  from  the  origin  to  a  given  point  of  stress- strain  curve.  The  results 
are  shown  in  the  Table  3.1. 
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For  those  materials  having  a  nonlinear  elastic  region,  use  of  the  strain  offset  method  is 
not  possible,  and  the  usual  practice  is  to  define  the  yield  strength  as  the  stress  required  to 
produce  some  amount  of  strain  [14],  In  this  work  the  yield  strength  is  defined  as  the 
stress  required  to  produce  a  strain  equal  to  0.003. 

Ao  =  WT  (3.3) 

a’~i  (3'4) 

pf 

<*f  =~T  (3-5) 

A, 

E  =  ^~  (3.6) 

Ae 

Table  3.2  Tensile  test  results 


rr.  .  ,  T ,  .  .  S  tan  dard  Deviation  .. 

Coefficient  of  Variation  = -  (3.7) 

Average 
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Table  3.  3  Results  average,  standard  deviation,  and  coefficient  of  variation 


E 

(psi) 

s  y 

(psi) 

sf 

(psi) 

Elongation 

(%) 

Average 

5.720  x  106 

12946 

16775 

0.475 

Standard 

Deviation 

7.330  x  105 

755 

3482 

0.108 

Coefficient  of 
Variation 

0.13 

0.06 

0.20 

0.23 
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Figure  3.2  Stress- strain  curves  for  copper  plated  film. 
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3.2  Wetting  Angle 


The  wetting  or  contact  angle  between  a  fluid  and  solid  boundary  varies  for  a  given  fluid 
with  the  type  of  solid  material,  and  the  roughness  and  cleanliness  of  the  surface.  Methods 
typically  used  to  determine  the  wetting  angle  are  the  titling  plate  (Spreece,  1957),  the 
sessile  drop  (Bickerman,  1941),  the  porous  plug  method  (Bartell,  1934),  and  the  wetting 
balance  method  (Guastalla,  1957)  [3].  The  method  used  in  this  study  is  the  sessile  drop. 

3.2.1  Equipment  and  procedure 


Figure  3.3  Contact  Angle  Meter. 


The  equipment  used  is  a  Tantec  Contact  Angle  Meter  model  CAM- Micro  illustrated  in 
Figure  3.3  and  the  procedure  is  the  following  [15]: 
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1.  Focus  the  needle  and  deposit  a  10  units  drop  over  the  surface. 

2.  Using  the  vertical  adjustment  knob,  position  the  droplet’s  image  so  its  base  is  on  the 
horizontal  cross- line. 

3.  Slide  the  measuring  screen  horizontally  until  the  left  edge  of  the  droplet’s  image 
touches  the  vertical  cross- line. 

4.  Determine  the  location  of  the  droplet’s  apex. 

5.  Rotate  the  clear  circular  piece  on  the  measuring  screen  until  the  hairline  crosses  the 
apex. 

6.  Read  the  contact  angle  on  the  scale. 


Table  3.4  Contact  angle  measurement. 


Surface 

Liquid 

Surface  Conditions 

Contact  Angle 
(Degrees) 

Plated  copper 

Tap  water 

Before  copper  cleaning 

90 

Plated  copper 

Tap  water 

After  copper  cleaning 

74 

Plated  copper 

DI  water 

Before  copper  cleaning 

75 

Plated  copper 

DI  water 

After  copper  cleaning 

56 

PWB  copper  layer 

Tap  water 

Before  copper  cleaning 

100 

PWB  copper  layer 

Tap  water 

After  copper  cleaning 

78 

PWB  copper  layer 

DI  water 

Before  copper  cleaning 

80 

PWB  copper  layer 

DI  water 

After  copper  cleaning 

60 

FR404  Epoxy 

Tap  water 

Clean 

120 

FR404  Epoxy 

DI  water 

Clean 

104 
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3.3  Heat  pipe  charging 


The  heat  pipe  was  charged  using  the  filling  method  based  on  a  micro- syringe  and  T- 
junction  adapter  developed  by  Y.  Cao,  M.  Gao,  and  E.  Pinilla  [16].  Prior  to  charging  the 
heat  pipe,  it  must  be  verified  that  the  heat  pipe  structure  is  hermetic. 

3.3.1  Leak  detection 

Several  methods  are  available  which  can  detect  leaks  of  different  magnitudes.  Very  large 
leaks  can  be  found  by  pressurizing  the  heat  pipe  with  compressed  gas  and  coating  the 
surface  with  a  soap  solution,  which  bubble  at  the  site  of  the  leak.  Smaller  leaks  are  notice 
when  the  pipe  cannot  be  pumped  down  sufficient  with  a  mcuum- roughing  pump.  Still 
smaller  leaks  will  cause  the  pressure  within  the  high-vacuum  pump  to  not  reach  its 
ultimate  capacity.  Therefore,  a  series  of  procedures  is  needed  to  determine  the  magnitude 
of  the  leak,  and  then  the  location  of  the  hole  in  the  container.  The  magnitude  of  the  leak 
will  be  quantified  by  the  pressure  seen  in  the  leaking  pipe  if  attached  to  a  vacuum  system, 
instead  of  the  standard  volume  of  air  per  second  (std  cm3/s).  This  is  done  since  most 
experiments  have  access  to  vacuum  equipment,  but  do  not  have  a  leak  rate  measurement 
system.  With  this  method,  the  pressure  inside  containers  with  gross  leaks  ranges  from 
atmospheric  to  1  torr,  low  vacuum  leaks  range  from  1  torr  to  1 0‘3  torr,  and  high- vacuum 
leaks  from  10‘3  to  10'9  torr.  High  vacuum  leaks  can  be  detected  with  a  spectrometer  leak 
detector. 
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3.3.2  Filling  the  Heat  Pipe  with  Working  Fluid 


After  passing  the  leakage  testing  the  laminated  heat  pipe  was  charged  by  using  syringe 
method.  In  this  method,  a  unique  micro- syringe,  designed  primarily  for  precision 
measurement  and  displacement  of  liquid  or  gaseous  samples  into  chromatographic 
analyzers  was  used.  The  syringe  consists  of  a  turn  valve  with  a  Teflon  seal,  a  leak-tight 
Teflon  plunger  tip,  and  a  removable  needle.  This  type  of  syringe  is  recommended  for  use 
in  either  a  high  vacuum  or  with  pressures  up  to  250  psi,  with  gases  or  liquids.  The  turn 
valve  can  be  used  to  isolate  the  working  fluid  stored  in  the  barrel  from  the  vacuum  during 
evacuating  of  the  heat  pipe.  The  plunger  makes  flush  contact  with  the  bottom  of  the 
syringe  barrel  when  fully  depressed,  virtually  eliminating  any  dead  volume. 

The  filling  setup  is  schematically  shown  in  Figure  3.  4. 


1 - Microsyringe  4-Heat  pipe 

2- Teflon  seal  5-To  vacuum  system 

3- Polyurethane  tubing 

Figure  3.  4  Micro- syringe  filling  setups 


The  heat  pipe  is  connected  to  the  syringe  using  the  polyurethane/Teflon  adapter  through 
which  the  needle  is  inserted,  it  is  completed  through  a  T-junction.  The  filling  procedure 
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for  the  system  incorporating  the  micro- syringe  and  the  T-junction  can  be  outlined  as 
follows: 

(1)  Fill  the  syringe  with  the  working  liquid. 

(2)  Lock  in  the  liquid.  Replace  the  needle  with  a  dry  one  and  lock  in  the  working 
liquid  in  the  barrel  by  screwing  the  nosepiece  with  at  least  two  turns  left  to 
tighten.  At  this  time,  there  is  no  liquid  in  the  needle,  and  the  liquid  is  isolated 
from  the  ambient. 

(3)  Evacuate  the  heat  pipe.  Connect  the  syringe  with  the  heat  pipe  through  the  T 
junction  and  connect  the  T-junction  to  the  vacuum  system,  until  the  desired 
vacuum  is  reached.  The  vacuum  inside  the  heat  pipe  is  difficult  to  measure  and  is 
certainly  lower  than  the  system  vacuum  because  of  the  small  diameter  of  the 
filling  pipe.  Adequate  time  should  be  allowed  for  the  heat  pipe  to  reach  the 
desired  vacuum.  Usually,  3  x  10'5  torr  of  vacuum  can  be  obtained. 

(4)  Fill  the  heat  pipe.  Screw  the  needle  nosepiece  all  the  way  and  depress  the  plunger 
down  to  the  dead  point  to  displace  the  working  liquid  into  the  heat  pipe  via  the 
filling  tube. ).  The  laminated  heat  pipe  was  charged  of  0.1  ml  of  distilled  water. 

(5)  Finally,  pinch  and  seal  the  heat  pipe  filling  tube. 

3.4  Heat  Pipe  Performance 

Steady  state  testing  and  evaluation  procedures  for  heat  pipes  and  thermosyphons  are 
typically  designed  to  measure  the  maximum  heat  transport  capacity,  the  end-to-end 
temperature  drop,  the  thermal  resistance,  and  /  or  the  effect  of  variations  in  the  local 
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gravitational  acceleration  on  the  overall  thermal  performance.  For  low  or  moderate 
temperature  heat  pipes  of  the  type  typically  encountered  in  the  thermal  control  of 
electronics,  all  of  these  measurements  can  be  accomplished  using  a  single,  relatively 
simple  test  facility  and  straightforward  test  procedures.  Figure  3.5  illustrates  a  steady 
state  test  facility  of  the  type  that  would  normally  be  used  for  low  or  moderate  temperature 
heat  pipes.  As  shown,  the  test  facility  consists  of  four  principal  subsystems:  the  test  pipe, 
a  heat  source,  a  heat  sink,  and  a  temperature  measurement  system. 


Input  Power  =  VI 


Electrical  Heating  Printed  Wiring  Board 

Figure  3.5  Heat  pipe  thermal  testing  setting 


3.4.1  Test  Facility 

Because  of  the  simplicity  and  accuracy  with  which  the  thermal  power  applied  to  the  test 
pipe  can  be  measured,  electrical  resistance  heating  is  the  most  common  type  of  heating 
system  used.  Other  methods,  including  induction  heating,  heated  constant  temperature 
baths,  and  radiation,  have  also  been  used  successfully  [1].  Regardless  of  the  method, 
what  is  needed  is  a  means  by  which  the  temperature  and  power  input  into  a  specific 
region  of  the  heat  pipe  can  be  carefully  monitored  and  controlled.  For  the  case  of 
electrical  resistance  heating,  the  electrical  energy  input  to  the  heat  pipe  or  joule  heating 
can  be  measured  either  by  measuring  the  input  voltage  and  current  or  by  measuring  the 
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voltage  and  using  the  known  electrical  resistance  of  the  heating  element.  The  electrical 
power  input  to  the  system  can  be  controlled  by  means  of  a  simple  variable  power 
transformer. 

In  addition  to  a  means  for  adding  power  to  the  evaporator  section  of  the  heat  pipe,  a 
system  is  also  required  by  which  the  heat  can  be  rejected.  Again,  there  are  several 
alternatives  available,  including  radiation,  latent  heat  rejection  to  a  spray  cooling  system, 
rejection  to  a  liquid- cooled  heat  exchanger,  and  free  convection  to  the  surrounding 
atmosphere  [1].  The  latter  of  these  was  used  in  this  study  because  represents  the  most 
practical  alternative. 

To  determine  the  transport  characteristics,  it  is  also  necessary  to  evaluate  the  axial 
temperature  distributions  along  the  length  of  the  heat  pipe.  These  measurements  were 
made  on  five  different  points  along  the  longitudinal  axis  on  the  external  case  of 
embedded  heat  pipe  (see  Figure  3.6).  Thermocouples  in  conjunction  with  a  computerized 
data  acquisition  system  were  used  to  monitor  the  power  input,  heat  rejected,  and 
temperature  distribution  throughout  the  tests.  While  the  temperature  gradient  along  the 
length  of  the  heat  pipe  may  be  quiet  small,  is  important  to  note  that  in  most  heat  pipe  test 
procedures,  the  absolute  temperature  measurements  are  not  as  critical  as  the  relative 
temperatures  differences.  For  this  reason,  the  use  of  thermocouples  with  a  relatively 
large  uncertainly  (±0.5  °F)  when  compared  with  the  temperature  gradients  being 
measured,  can  yield  results  that  are  quite  acceptable.  However,  care  should  be  taken  to 
ensure  that  at  any  given  temperature  all  of  the  thermocouples  would  yield  the  same  value. 
This  was  accomplished  by  verifying  the  thermocouples  and  data  acquisition  system  at  the 
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freezing  point  and  atmospheric  boiling  point  of  water  to  ensure  that,  for  a  given 
condition,  they  yield  similar  results  prior  the  installation  on  or  in  the  test  pipe.  Previous 
tests  have  shown  that  temperature  differences  can  be  measured  to  within  an  accuracy  of 
+0.05  °F  using  this  technique  [1], 


Figure  3.6  Temperature  points  distributions  along  the  length  of  the  heat  pipe. 


3.4.2  Test  Procedure 

Test  procedures  consisted  in  incrementing  the  input  power  to  the  evaporator  in 
predetermined  increments  and  determined  the  maximum  power  input  that  the  embedded 
heat  pipe  in  the  PWB  could  accept  without  failing.  The  increments  in  power  input  were: 
2.5  W,  5.6  W,  and  10  W.  The  same  tests  were  conducted  to  the  PWB  with  uncharged 
heat  pipe.  Figure  3.7  to  Figure  3.9  compare  the  temperature  distribution  along  heat  pipe 
length  for  charged  and  uncharged  heat  pipes  for  2.5  W,  5.6  W,  and  10  W  respectively. 
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Temperature  (F) 


■2.5  W-Charged  HP 
•2.5  W- Uncharged  HP 


Figure  3.7  Experimental  temperature  distribution  along  heat  pipe  length  with  2.5  W  of  input  power. 
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—*—5.6  W-Charged  HP 
5.6  W-Uncharged  HP 


Figure  3.8  Experimental  temperature  distribution  along  heat  pipe  length  with  5.6  W  of  input  power. 
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10  W-Charged  HP 
10  W-Uncharged  HP 


Figure  3.9  Experimental  temperature  distribution  along  heat  pipe  length  with  10  W  of  input  power. 
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Section  4 


Finite  Element  Modeling 

In  parallel  with  the  test,  a  finite  element  analysis  was  also  made  in  the  present  study  to 
determine  the  stresses  induced  in  the  substrate  by  thermal  gradients.  The  intent  was  to 
reproduce  the  thermal  gradients  found  in  the  experimental  thermal  performance  fitting 
the  heat  pipe  conductivity  in  the  finite  element  model. 

The  physical  structure  of  the  embedded  heat  pipe  was  modeled  for  both  heat  pipe 
operation  and  system  level  exposure  to  environmental  stresses,  such  as  temperature, 
using  a  FEA  software  (ABAQUS/CAE  6.1  from  HKS,  Inc).  The  material  properties  of 
the  laminated  material,  and  the  plated  metallization  were  determined  and  used  to  model 
the  mechanical  performance  (i.e.  thermally  induced  stresses  and  strains)  over  the 
environmental  stress  exposure.  For  this  purpose  the  following  two  FEA  analyses  were 
conducted: 

1.  A  3D  heat  transfer  analysis  to  determine  the  temperature  distribution  through  the 
whole  PWB  model. 

2.  A  2D  sequentially  coupled  thermal- stress  analysis  of  the  cross  section  at  a  selected 
point  of  the  embedded  heat  pipe. 

Because  ABAQUS  has  no  built-in  system  a  consistent  system  of  units  has  to  be  selected 
for  all  input  data  before  starting  to  define  the  models  [17].  The  system  selected  was  the 
US  Unit  (inch)  system. 

Length:  in 
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Force:  lbf 


Mass:  lbf-  s2/in 
Time:  s 
Stress:  lbf/in2 
Energy:  in- lbf 
Power:  in-lbf/s 
Density:  lbf  s2/in4 

4.1  3D  Heat  Transfer  Analysis  of  the  PWB  Model 

An  uncoupled  3D  heat  transfer  analysis  is  used  to  model  the  PWB  solid  body  with 
general,  constant  thermal  conductivity  irrespective  of  temperatures,  internal  energy 
generated  in  the  heater,  and  quite  general  convection  and  radiation  boundary  conditions. 

Two  different  analyses  were  performed  in  this  model.  The  first  one  simulating  a  charged 
heat  pipe  and  the  second  one  simulating  an  uncharged  heat  pipe.  The  only  difference 
between  these  two  analyses  is  the  thermal  conductivity  of  the  vapor  space  part  as  shown 
in  Table  4.1 

4.1.1  Assumptions 

The  numerical  analysis  was  conduced  under  the  following  assumptions: 

(1)  Heat  transfer  coefficient  from  test  printed  circuit  board  to  air  has  a  constant  value 
irrespective  of  location  and/or  shape. 

(2)  The  ambient  air  temperature  is  constant. 
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(3)  A  printed  circuit  board  model  has  a  shape  and  dimensions  shown  in  Figure  2.  5  and 
Figure  2.6. 

(4)  Every  member  has  a  constant  thermal  conductivity  irrespective  of  temperatures. 

(5)  Contact  resistance  is  constant  between  the  heating  element  and  the  printed  circuit 
board. 

(6)  The  axial  temperature  testing  points  along  the  length  of  the  heat  pipe  were  located  at 
the  same  positions  as  the  experimental  testing  points  (see  Figure  3.6). 

4.1.2  Model  Definition 

The  defining  parameters  of  the  model  include  shape  and  geometry  information  as  well  as 
material  property  data.  All  of  these  components  are  modeled  using  eight-node 
convection  diffusion  brick  linear  element.  The  typical  material  properties  were  used, 
except  for  the  conductivity  of  PWB  material.  Since  the  typical  conductivity  of  the  FR404 
epoxy  is  not  the  same  as  the  PWB  part  after  lamination  due  the  copper  layer  in  the  top 
and  bottom  surfaces  of  each  laminate  layer.  The  conductivity  of  the  PWB  material  was 
fit  to  meet  the  experimental  temperature  distribution. 

Table  4.1  lists  the  material  properties  of  the  different  parts  of  the  heat  pipe  assembly 
including  the  two  different  sets  of  materials  properties  for  the  Vapor  Space  part,  one  for 
the  charged  heat  pipe  analysis  and  other  one  for  the  uncharged  heat  pipe  analysis.  These 
parts  can  be  seen  in  Figure  4.1. 


51 


Vapor  Space 


PWB 


Copper  Layers 


Heater 


Figure  4.1  Printing  Wiring  Board  model. 


Table  4.1  3D  Thermal  Model  material  properties. 


Conductivity 

[W/m°K] 

Density 

[g/cm3] 

Specific  Heat 
[J/kg°K] 

Copper  Layer 

398 

8.96 

385 

PWB 

X 

7.61 

Y 

0.48 

z 

7.61 

1.90 

1150 

Vapor  Space 
(Charged  Heat  Pipe 
Analysis) 

X 

80 

Y 

400 

Z 

80 

8.96 

385 

Vapor  Space 
(Uncharged  Heat  Pipe 
Analysis) 

X 

0 

Y 

0 

Z 

0 

8.96 

385 

Heater 

398 

8.96 

385 
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4.1.3  Element  Selection 


The  element  selected  for  the  all  mesh  model  components  is  the  linear  hexahedron 
element,  8-node  convection  diffusion  brick,  ABAQUS  type  DCC3D8.  This  element 
holds  temperature  as  active  degree  of  freedom  [17]. 

4. 1.3.1  Surface  Contact  Interaction 

The  Surface  Contact  Interaction  in  this  model  defines  the  thermal  interactions  between 
regions  of  the  model  and  between  a  region  of  the  model  and  its  surroundings. 

The  Tie  Surface  Interaction  was  used  to  fuse  together  all  part-instance  faces  in  contact 
even  though  the  meshes  created  on  the  surfaces  may  be  dissimilar.  Creating  a  tie  surface 
interaction  is  analogous  to  including  the  *CONTACT  PAIR,  TIED  option  h  a  solver 
input  file. 

The  Thermal  Film  Condition  was  used  to  define  the  cooling  due  to  convection  of  all 
surfaces  of  the  model  expose  to  the  surrounding.  A  film  condition  interaction  is 
analogous  to  including  the  *SFILM  option  in  a  solver  input  file. 

Radiation  to  Ambient  Interaction  used  in  this  model  describes  heat  transfer  between  the 
nonconcave  surface  of  the  model  and  a  nonreflecting  environment  due  to  radiation.  A 
radiation  interaction  is  analogous  to  including  the  *SRADIATE  option  in  a  solver  input 
file. 


4.1.4  Mesh  Discretization 

ABAQUS  /CAE  can  use  a  variety  of  meshing  techniques  to  mesh  models  of  different 
topologies:  structured  meshing,  swept  meshing,  free  meshing.  The  different  meshing 
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techniques  provide  varying  levels  of  automation  and  user  control.  Figure  4.2  displays  the 


meshed  model. 


Figure  4.2  Meshed  Model. 

Structure  meshing  was  used  since  it  gives  most  control  over  the  mesh  because  it  applies 
pre-established  mesh  patterns  to  particular  model  topologies.  The  PWB- 1  part- instance 
was  too  complex  to  be  meshed  using  pre-established  mesh  patterns  directly. 
Consequently,  it  had  to  be  partitioned  into  simple  regions  with  topologies  for  which 
structured  meshing  patterns  exits.  The  total  number  of  elements  is  3435  and  nodes  7276. 
Table  4.2  shows  the  mesh  break  down  by  part- instances. 
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Table  4.2  Part- instance  mesh  break  down. 


PART- 

INSTANCE 

ELEMENT 

TYPE 

NODES 

ELEMENTS 

Copper  Layer-1 

Linear  Hexahedron 
DCC3D8 

2806 

1320 

Copper  Layer-2 

Linear  Hexahedron 
DCC3D8 

2880 

1357 

PWB-1 

Linear  Hexahedron 
DCC3D8 

384 

214 

Vapor  Space -1 

Linear  Hexahedron 
DCC3D8 

198 

84 

Heater-1 

Linear  Hexahedron 
DCC3D8 

1008 

460 

TOTAL 

7276 

3435 

4.1.5  Loading  Conditions 

The  loading  condition  attempt  to  simulate  the  heat  generated  in  the  heater.  Three  heat 
inputs  levels  will  be  analyzed:  2.5  W  (22.121  in  lbf/s),  5.6  W  (49.551  in  lbf/s),  and  10  W 
(88.485  in  lbf/s). 

The  *DFLUX  option  was  used  to  specify  the  body  fluxes  (flux  per  unit  volume)  in  the 
heater.  Since  the  heater  volume  is  0.06816  in3  the  corresponding  flux  per  unit  volume  for 
the  three  different  heat  inputs  are:  36.678  W/in3  (324.547  in  lbf/s  in3),  82.160  W/in3 
(726.985  in  lbtfs  in3),  and  146.714  W/in3  (1298.188  in  lbf/s  in3). 

4.1.6  Type  of  Analysis 

Uncoupled  heat  transfer  analysis  is  a  pure  heat  transfer  analysis  using  heat  transfer 
elements  and  the  *HEAT  TRANFER  procedure.  This  is  linear,  transient,  heat  transfer 
analysis  used  to  model  solid  body  heat  conduction  with  general,  constant  thermal 
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conductivity  irrespective  of  temperatures,  internal  energy  generated  in  the  heater,  and 
quite  general  convection  and  radiation  boundary  conditions. 
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Temperature  (F) 


2.5  W-Charged  HP 
2.5  W- Uncharged  HP 


Figure  4.3  FEM  temperature  distribution  along  heat  pipe  length  with  2.5  W  of  input  power. 
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Figure  4.4  Temperature  distribution  top  view  for  2.5  W  of  input  power:  (a)  charged  heat 
pipe,  (b)  uncharged  heat  pipe. 
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Figure  4.5  Temperature  distribution  bottom  view  for  2.5  W  of  input  power:  (a)  charged 
heat  pipe,  (b)  uncharged  heat  pipe. 
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—*—5.6  W-Charged  HP 
5.6  W-Uncharged  HP 


Figure  4.6  FEM  temperature  distribution  along  heat  pipe  length  with  5.6  W  of  input  power. 
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Figure  4.7  Temperature  distribution  top  view  for  5.6  W  of  input  power:  (a)  charged  heat 
pipe,  (b)  uncharged  heat  pipe. 
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Figure  4.8  Temperature  distribution  bottom  view  for  5.6  W  of  input  power:  (a)  charged 
heat  pipe,  (b)  uncharged  heat  pipe. 
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10  W-Charged  HP 
10  W-Uncharged  HP 


Figure  4.9  FEM  temperature  distribution  along  heat  pipe  length  with  10  W  of  input  power. 
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Figure  4.  10  Temperature  distribution  top  view  for  10  W  of  input  power:  (a)  charged 
heat  pipe,  (b)  uncharged  heat  pipe. 
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Figure  4.  1 1  Temperature  distribution  bottom  view  for  10  W  of  input  power:  (a)  charged 
heat  pipe,  (b)  uncharged  heat  pipe. 
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Figure  4.  12  Experimental  and  FEM  temperature  distributions  along  charged  heat  pipe  with  2.5  W  of  input  power. 
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Experimental 


—“—FEM 


Figure  4.13  Experimental  and  FEM  temperature  distributions  along  charged  heat  pipe  with  5.6  W  of  input  power. 
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Figure  4.14  Experimental  and  FEM  temperature  distributions  along  charged  heat  pipe  with  10  W  of  input  power. 
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Figure  4.15  Experimental  and  FEM  temperature  distributions  along  uncharged  heat  pipe  with  2.5  W  of  input  power. 
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Figure  4.16  Experimental  and  FEM  temperature  distributions  along  uncharged  heat  pipe  with  5.6  W  of  input  power. 
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Figure  4.17  Experimental  and  FEM  temperature  distributions  along  uncharged  heat  pipe  with  10  W  of  input  power 
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4.2  2D  Sequentially  Coupled  Thermal -Stress  Analysis  of  the  PWB  Cross  Section 


A  sequentially  coupled  heat  transfer  analysis  is  used  when  the  stress/deformation  field  in 
a  structure  depends  on  the  temperature  field  in  that  structure,  but  the  temperature  field 
can  be  found  without  knowledge  of  the  stress/deformation  response.  A  sequentially 
coupled  heat  transfer  analysis  is  performed  by  first  conducting  a  pure  heat  transfer 
analysis  and  then  a  stress/deformation  analysis. 

Nodal  temperature  are  stored  as  a  function  of  time  in  the  heat  transfer  results  (.fil)  file  by 
using  the  *NODE  FILE  option  and  specifying  output  variable  NT. 

To  define  the  temperature  field  at  different  times  in  the  stress  analysis  ^TEMPERATURE 
option  is  used  within  the  stress  analysis  step  definition  to  read  temperatures  at  the  nodes 
stored  in  the  heat  transfer  results  file.  This  interface  assumes  that  node  numbers  are  the 
same  for  corresponding  nodes  in  the  stress  analysis  mesh  and  in  the  heat  transfer  analysis 
mesh. 

The  Vapor  Space- 1  and  Heater- 1  part- instances  were  omitted  for  stress/displacement 
analysis  to  avoid  their  interaction  with  the  rest  of  the  structural  elements  of  the  model. 
Nodes  from  Vapor  Space- 1  and  Heater- 1  part- instances  were  removed  for  the  stress 
problem  using  the  *MODEL  CHANGE  option.  When  the  results  file  is  read  during  the 
stress  analysis,  temperature  at  nodes  that  are  not  in  the  mesh  are  ignored. 
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4.2.1  Assumptions 


The  numerical  analysis  was  conduced  under  the  following  assumptions: 

(1)  The  stress/deformation  field  depends  on  the  temperature  field,  but  the  temperature 
field  can  be  found  without  knowledge  of  the  stress/deformation  response. 

(2)  Inertia  effects  are  neglected. 

(3)  Ignores  time- dependent  material  effects  (creep,  swelling,  viscoelasticity). 

(4)  Heat  transfer  coefficient  from  test  printed  circuit  board  to  air  has  a  constant  value 
irrespective  of  location  and/or  shape. 

(5)  The  ambient  air  temperature  is  constant. 

(6)  A  printed  circuit  board  cross  section  model  has  a  shape  dimensions  as  shown  in 
Figure  2.6. 

(7)  Every  member  has  a  constant  thermal  conductivity  irrespective  of  temperatures. 

(8)  Contact  resistance  is  constant  between  the  heating  element  and  the  printed  circuit 
board. 

4.2.2  Model  Definition 

In  the  heat  transfer  analysis  all  components  are  modeled  using  a  four- node  convection 
diffusion  quadrilateral  linear  element,  while  in  the  stress/deformation  analysis  all  the 
components  are  modeled  using  a  four-node  bilinear  plane  strain  quadrilateral  linear 
element.  Table  4.3  lists  the  material  properties  of  the  different  parts  of  the  heat  pipe 
assembly.  These  parts  can  be  seen  in  Figure  4.  18. 
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Table  4.3  2D  Thermal- Stress  Model  material  properties. 


Conductivity 

[W/m°K] 

Density 

[g/cnrj 

Specific  Heat 
[J/kg°K] 

Modulus  of 
Elasticity 
GPa 

Poisson’s 

Ratio 

Coefficient  of 
Thermal  Expansion 
1<F  (°C)-' 

Copper 

Layer 

398 

8.96 

385 

110 

0.35 

16.5 

FR-4 

X 

7.61 

Y 

0.48 

Z 

7.61 

1.90 

1150 

17.4 

0.19 

X 

14 

Y 

175 

Z 

13 

Vapor 

Space 

X 

80 

Y 

400 

Z 

80 

8.96 

385 

110 

0.35 

16.5 

Heater 

398 

8.96 

385 

110 

0.35 

16.5 

Copper 

Plating 

398 

8.96 

385 

39.36 

0.35 

16.5 

Figure  4.  18  PWB  Cross  Section  Model. 

4.2.3  Element  Selection 

The  element  selected  in  the  heat  transfer  analysis  for  the  all  mesh  model  components  is 
the  4-node  convection  diffusion  quadrilateral  linear  element  ABAQUS  type  DCC2D4. 
This  element  holds  temperature  as  active  degree  of  freedom  [17]. 
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The  element  selected  for  the  all  mesh  model  components  in  the  stress/displacement 
analysis  is  the  4-node  bilinear  plane  strain  quadrilateral  linear  element  ABAQUS  type 
CPE4R.  This  element  holds  1,  2  (Ux,  Uy)  as  active  degree  of  freedom  [17]. 

4.2.3. 1  Surface  Contact  Interaction 

The  Tie  Surface  Interaction  was  used  to  fuse  together  all  part- instances  surfaces  in 
contact  even  though  the  meshes  created  on  the  surfaces  may  be  dissimilar.  Creating  a  tie 
surface  interaction  is  analogous  to  including  the  *CONTACT  PAIR,  TIED  option  in  a 
solver  input  file. 

In  the  heat  transfer  analysis  the  Thermal  Film  Condition  is  used  to  define  the  cooling  due 
to  convection  of  all  surfaces  of  the  model  expose  to  the  surrounding.  A  film  condition 
interaction  is  analogous  to  including  the  *SFILM  option  in  a  solver  input  file. 

Radiation  to  Ambient  Interaction  is  used  in  the  heat  transfer  analysis  to  model  heat 
transfer  between  the  nonconcave  surface  of  the  model  and  a  nonreflecting  environment 
due  to  radiation.  A  radiation  interaction  is  analogous  to  including  the  *SRADIATE 
option  in  a  solver  input  file. 

4.2.4  Mesh  Discretization 

Structure  meshing  was  used  since  it  gives  most  control  over  the  mesh  because  it  applies 
pre-established  mesh  patterns  to  particular  model  topologies.  The  FR-4-1  part- instance 
was  too  complex  to  be  meshed  using  pre-established  mesh  patterns  directly. 
Consequently,  it  had  to  be  partitioned  into  simple  regions  with  topologies  for  which 
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structured  meshing  pattern  exits.  The  total  number  of  elements  is  1550  and  nodes  2443. 
Table  4.4  shows  the  mesh  break  down  by  part- instances. 


i 


Figure  4.19  Meshed  Model 


Table  4.4  Part-  instance  mesh  break  down 


PART-INSTANCE 

HEAT  TRANSFER 
ELEMENT  TYPE 

STRESS/DISPLACEMENT 
ELEMENT  TYPE 

NODES 

ELEMENTS 

Copper  Layer-1 

Linear  quadrilateral  DCC2D4 

Linear  quadrilateral  CPE4R 

452 

225 

Copper  Layer-2 

Linear  quadrilateral  DCC2D4 

Linear  quadrilateral  CPE4R 

452 

225 

FR-4-1 

Linear  quadrilateral  DCC2D4 

Linear  quadrilateral  CPE4R 

746 

637 

Vapor  Space -1 

Linear  quadrilateral  DCC2D4 

Linear  quadrilateral  CPE4R 

180 

144 

Heater-1 

Linear  quadrilateral  DCC2D4 

Linear  quadrilateral  CPE4R 

87 

56 

Copper  Plating 

Linear  quadrilateral  DCC2D4 

Linear  quadrilateral  CPE4R 

526 

263 

TOTAL 

2443 

1550 

76 


4.2.5  Loading  Conditions 


The  loading  condition  attempt  to  simulate  the  heat  generated  in  the  heater.  Three  heat 
inputs  levels  will  be  analyzed:  2.5  W  (22.121  in  lb./s),  5.6  W  (49.551  in  lb./s),  and  10  W 
(88.485  in  lb./s).  The  *DSFLUX  option  was  used  to  specify  the  distributed  surface 
fluxes  in  the  heater. 

4.2.6  Initial  Conditions 

By  default  the  initial  temperature  of  all  nodes  is  zero.  The  TYPE=TEMPERATURE 
parameter  is  used  on  the  *INITIAL  CONDITIONS  option  to  specify  nonzero  initial 
temperatures. 

4.2.7  Boundary  Conditions 

Only  zero- valued  boundary  conditions  can  be  prescribed  as  model  data.  The  data  lines 
can  be  specified  using  either  “type”  or  “direct”  format.  The  “Type”  format  is  a  way  of 
conveniently  specifying  common  types  of  boundary  conditions  in  stress/displacement 
analyses. 

For  “type”  format  the  first  entry  on  the  data  line  is  the  node  set  or  node  number  to  which 
the  boundary  conditions  apply.  The  following  entries  specify  the  degrees  of  freedom  to 
be  restrained  at  that  node  or  node  set.  Either  a  single  degree  of  freedom  or  the  first  and 
last  of  a  range  of  degree  of  freedom  can  be  specified  directly  on  the  data  lines  of  the 
*BOUNDARY  option.  The  bottom  comers  of  the  FR-4  part  instance  are  constrained  in 
d.o.f  U2  (vertical  displacement)  (see  Figure  4.20). 
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Figure  4.20  2D  Stress  Analysis  Model  Boundary  Conditions. 


4.2.8  Type  of  Analysis 

The  sequentially  coupled  heat  transfer  analysis  is  performed  running  first  an  uncoupled 
heat  transfer  analysis  followed  by  a  stress/deformation  analysis  run.  Nodal  temperature 
is  stored  at  the  heat  transfer  analysis  in  the  heat  transfer  results  (.fil)  file  and  subsequently 
it  is  read  at  the  stress/displacement  analysis. 

Uncoupled  heat  transfer  analysis  is  a  pure  heat  transfer  analysis  using  heat  transfer 
elements  and  the  *HEAT  TRANFER  procedure.  This  is  linear,  steady  state,  heat  transfer 
analysis  used  to  model  solid  body  heat  conduction  with  general,  constant  thermal 
conductivity  irrespective  of  temperatures,  internal  energy  generated  in  the  heater,  and 
quite  general  convection  and  radiation  boundary  conditions. 
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Stress/deformation  analysis  is  linear,  steady  state,  static  analysis  and  the  *STATIC 
procedure  is  used  to  perform  this  analysis.  The  response  is  caused  by  temperatures 
varying  with  time  based  on  the  previous  transient  heat  transfer  run. 
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Figure  4.21  Temperature  Distribution  at  point  T1  with  2.5  W  of  input  power. 
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Figure  4.22  Mises  stress  at  point  Tl  with  2.5  W  of  input  power. 
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Figure  4.23  S22  stress  at  point  T1  with  2.5  W  of  input  power. 
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Figure  4.24  Temperature  Distribution  at  point  T3  with  2.5  W  of  input  power. 
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Figure  4.25  Mises  stress  at  point  T3  with  2.5  W  of  input  power. 
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Figure  4.26  S22  stress  at  point  T3  with  2.5  W  of  input  power. 
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Figure  4.  27  Temperature  Distribution  at  point  T5  with  2.5  W  of  input  power. 
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Figure  4.  28  Mises  stress  at  point  T5  with  2.5  W  of  input  power. 
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Figure  4.29  S22  stress  at  point  T5  with  2.5  W  of  input  power. 
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Figure  4.30  Temperature  Distribution  at  point  T1  with  5.6  W  of  input  power. 


89 


V™por1 1  O  DR:  r.:  ktraq  EHt_f ite  UPC  irasSetf_S1  ress_6_GW  jrwiti 


Sf  Hises 
(Ave.  Crm. : 

a+Z  _  5-63  e  #03 
+£.<i74*f03 
+2.38-ieK)3 
+Z.Z9*e#03 
+2.204**03 
+2  _  lHe*-03 
+2 . Q££et03 
+1.935**03 
+1 . 345 e *03 
+1 . ?55et03 
— h  +1  -  S  65  ft  *0  3 
+1 .  1 76  s? +03 
+1.486e*03 


£ 


1 


Figure  4.31  Mises  stress  at  point  T1  with  5.6  W  of  input  power. 
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Figure  4.32  S22  stress  at  point  T1  with  5.6  W  of  input  power. 
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Figure  4.33  Temperature  Distribution  at  point  T3  with  5.6  W  of  input  power. 
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Figure  4.34  Mises  stress  at  point  T3  with  5.6  W  of  input  power. 
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Figure  4.35  S22  stress  at  point  T3  with  5.6  W  of  input  power. 
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Figure  4.  36  Temperature  Distribution  at  point  T5  with  5.6  W  of  input  power. 
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Figure  4.37  Mises  stress  at  point  T5  with  5.6  W  of  input  power. 
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Figure  4.  38  S22  stress  at  point  T5  with  5.6  W  of  input  power. 
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Figure  4.39  Temperature  Distribution  at  point  T1  with  10  W  of  input  power. 
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Figure  4.40  Mises  stress  at  point  T1  with  10  W  of  input  power. 
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Figure  4.41  S22  stress  at  point  T1  with  10  W  of  input  power. 
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Figure  4.42  Temperature  Distribution  at  point  T3  with  10  W  of  input  power 
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Figure  4.43  Mises  stress  at  point  T3  with  10  W  of  input  power. 
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Figure  4.44  S22  stress  at  point  T3  with  10  W  of  input  power 


103 


Vpwiprarl  e  1  O  BB:  r-'-  tebnq  us_f  ifes/HPC  ™sS*?,,rgl_riiwima  l_  i  0 W_  Pftml^jnrlh 


MTU 


+1.7SSet02 
+1_747«*0Z 
+l_?33e*0£ 
+1 . 730ef02 
+1_7ZZ«#0Z 
+1.7l-fe^0E 
■+1.706ef02 
+1  _  G-S*7e#0,Z 

+1  .€89*mz 
+1 .  681e402 
+1 _  67ZetOZ 

+I.e6-U*Q£ 

+1.6S6ef 02 


Figure  4.45  Temperature  Distribution  at  point  T5  with  10  W  of  input  power. 
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Figure  4.46  Mises  stress  Distribution  at  point  T5  with  10  W  of  input  power 
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Figure  4.47  S22  stress  Distribution  at  point  T5  with  10  W  of  input  power 


106 


Viewport;  1  ODR:  n:  tfhar}  ik  _f ilps  MPC  mssS*ipi_S1  nss_  1  flW  jw*i 


£,.  Hises 

(Ave.  Crm. :  7£M 


3 - 


Figure  4.48  Maximum  stress  in  the  copper  plating  with  10  W  of  input  power.. 
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Table  4.5  Maximum  temperature  and  stresses  in  the  substrate  (FR-4)  cross-section  with 
2.5  W  of  input  power. 


Cross  Section 

Temperature 

Mises  Stress 

Sll 

S22 

at  Point 

(°C) 

(psi) 

(psi) 

(psi) 

T1 

51 

1234 

319 

550 

T2 

44 

980 

253 

437 

T3 

41 

807 

208 

360 

Table  4.6  Maximum  temperature  and  stresses  in  the  substrate  (FR-4)  cross-section  with 
5.6  W  of  input  power. 


Cross  Section 

Temperature 

Mises  Stress 

Sll 

S22 

at  Point 

(°C) 

(psi) 

(psi) 

(psi) 

T1 

81 

2563 

666 

1142 

T2 

68 

1999 

518 

891 

T3 

59 

1623 

420 

724 

Table  4.  7  Maximum  temperature  and  stresses  in  the  substrate  (FR-4)  cross-section  with 
10  W  of  input  power. 


Cross  Section 

Temperature 

Mises  Stress 

Sll 

S22 

at  Point 

(°C) 

(psi) 

(psi) 

(psi) 

T1 

113 

3970 

1035 

1768 

T2 

94 

3137 

817 

1397 

T3 

79 

2520 

655 

1123 
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Figure  4.49  Substrate  maximum  Mises  stresses  along  embedded  heat  pipe  length. 
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Appendix  A 
PTH  Chemistry 

Wear  protective  clothing  and  face  shield  when  handling  chemicals.  Always  add  acid  to 
water,  slowly;  never  the  reverse.  Fuming  and  heat  will  occur  when  mixing.  Avoid 
breathing  fumes 

Step  #  1  -  Chemical  Cleaner:  CCH-3,  3.5G,  PTH  Cleaner 
Pour  71ml.  of  water  into  a  beaker;  add  29  ml.  of  concentrate. 

Step  #  3  -  Copper  Cleaner:  AP-24  G,  Ammonium  Persulphate 

Pour  850  ml.  of  hot  water  into  a  beaker;  add  205  gr.  of  AP  crystals.  Stir  until  dissolved. 
To  this  solution,  add  18  ml.,  Sulfuric  Acid,  H2SO4,  Reagent  Grade  95-98  %. 

Step  #  5  -  Acid  Tread:  H2SO4,  Sulfuric  Acid,  Reagent  Grade  95-98  % 

Pour  857  ml  of  water  into  a  beaker;  slowly  add  143  ml  of  Sulfuric  Acid. 

Step  #  7  -  Acid  Treat:  HCL,  Hydrochloric  Acid,  Reagent  Grade  36.5-38.0% 

Pour  714  ml.  of  water  into  a  beaker;  add  286  ml.  of  hydrochloric  acid. 

Step  #  8  -  Activator:  ACT-1G,  Activator 

Used  as  supplied.  Approximately  100  ml.  required. 

Step  #  10  -  Accelerator:  ACC-3.5G,  Accelerator 

Pour  714  ml.  of  deionized  water  into  a  beaker;  add  286  ml.  of  concentrate. 
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Step  #  13  -  Electroless  Copper:  ECU-10G,  Electroless  Copper,  supplied  as  Part  “A” 


and  Part  “B”  concentrates 

Pour  100  ml.  of  Part  “A”  into  a  beaker;  add  800  ml.  of  deionized  water;  add  100  ml.  of 
Part  “B”  (Mixture  is  10%  Part  “A”,  80  %  DI  water,  10%  Part  “B”;  DO  NOT  mix  parts 
“A”  and  “B”  directly  together). 

Step  #  15  -  Acid  Treat:  H2SO4,  Sulfuric  Acid,  Reagent  Grade  95-98  % 

Pour  929  ml.  of  water  into  a  beaker;  slowly  add  71ml.  of  Sulfuric  Acid. 

Step  #  16  -  Copper  Electroplate:  CU-5G,  Copper  Electrolyte 
Use  as  supplied,  approximately  1000  ml. 
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